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Using analytical calculations and computer simulations we consider both the lateral diffusion of a 
membrane protein and the fluctuation spectrum of the membrane in which the protein is embedded. 
The membrane protein interacts with the membrane shape through its spontaneous curvature and 
bending rigidity. The lateral motion of the protein may be viewed as diffusion in an effective 
potential, hence, the effective mobility is always reduced compared to the case of free diffusion. 
Using a rigorous path-integral approach we derive an analytical expression for the effective diffusion 
coefficient for small ratios of temperature and bending rigidity, which is the biologically relevant 
limit. Simulations show very good quantitative agreement with our analytical result. The analysis of 
the correlation functions contributing to the diffusion coefficient shows that the correlations between 
the stochastic force of the protein and the response in the membrane shape are responsible for the 
reduction. 

Our quantitative analysis of the membrane height correlation spectrum shows an influence of the 
protein-membrane interaction causing a distinctly altered wave-vector dependence compared to a 
free membrane. Furthermore, the time correlations exhibit the two relevant timescales of the system: 
that of membrane fluctuations and that of lateral protein diffusion with the latter typically much 
longer than the former. We argue that the analysis of the long-time decay of membrane height 
correlations can thus provide a new means to determine the effective diffusion coefficient of proteins 
in the membrane. 

PACS numbers: 87.16.D-, 87.16.A-, 87.15.Vv, 87.16.dj, 05.40.-a 



I. INTRODUCTION 

Biomembranes are ubiquitous in life, mainly provid- 
ing spatial compartmentalisation. However, a membrane 
should not be viewed as a mere barrier between differ- 
ent compartments, but serves as a place where a whole 
variety of functions may take place, like ion or protein 
transport, signal transduction, etc. pQ. These functions 
come about through proteins that move along the mem- 
brane. From a physical perspective the lateral diffusion 
of the proteins and the shape changes of the membrane 
caused upon insertion of proteins are among the most 
interesting issues of these systems. 

The recent progress in experimental techniques to mea- 
sure lateral diffusion coefficients, like fluorescence cor- 
relation spectroscopy [2], single particle tracking [3], or 
fluorescence recovery after photobleaching 0], has re- 
vealed that many of the functions performed by pro- 
teins are crucially influenced by the diffusive behavior 
of the proteins [5] . Apart from the obvious biological rel- 
evance lateral protein diffusion is also very challenging 
from a theoretical perspective: Compared to diffusion in 
the bulk there is a subtlety in the hydrodynamic equa- 
tions describing the mobility in a two dimensional fluid, 
since the solution of the two-dimensional Navier-Stokes 
equation diverges. In order to overcome this so-called 
Stokes' paradox [B], Saffman and Delbriick 7 considered 
the mobility of a very thin, rigid object in a narrow al- 
most two-dimensional fluid layer that is surrounded on 
both sides by a further liquid. This work has received 
a lot of attention since it is relevant for lateral protein 
diffusion. While some experiments support their result 
[8HTU]. more recent observations for proteins cannot be 



explained by their theory [TTJ Q2] . 

Another aspect that makes diffusion interesting, par- 
ticularly in membranes, is that the membrane itself is 
subject to thermal fluctuations; thus the shape of the 
membrane is also constantly changing. Methods to anal- 
yse shape fluctuations of membranes include off-specular 
x-ray scattering [T3] and video microscopy [T21[T5]- In the 
latter method the contour of a vesicle is detected from 
optical microscopy records taken at successive timesteps. 
The changes in the contour provide information on the 
fluctuation spectrum that is used to deduce effective 
bending rigidities or surface tensions. In a very recent 
study Rodriguez-Garcfa et al. [16] identified the influence 
of the bilayer nature of a membrane as theoretical calcu- 
lations [T7] have previously predicted. The influence of 
the density of inclusions embedded in a lipid membrane 
on the effective bending rigidity was studied by Vitkova 
et al. [18j . In this work the peptide alamethicin was used 
as the inclusion. Bassereau and co-workers have studied 
more complicated systems consisting of membranes with 
inserted proteins and have been able to study the altered 
fluctuation spectrum of a membrane upon activation of 
the inserted bacteriorhodopsin proteins [THJ HO] • 

Theoretically, a bare membrane is well described as a 
continuous two-dimensional sheet with a bending rigid- 
ity and an effective surface tension. This model has been 
very successful in explaining a whole variety of exper- 
imentally observed membrane morphologies [21] . Like- 
wise, membrane shape fluctuations are well captured by 
this simple model as shown in video microscopy exper- 
iments [21 [15] . The insertion of additional proteins in 
a membrane requires an extension of this simple con- 
tinuous model to include the local interacion of a protein 
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with the membrane. While the influence of thermal mem- 
brane fluctuations on the interaction between inclusions 
has previously been considered in several studies [23] , 
the influence on lateral diffusion or the altered membrane 
height correlations is much less studied. In previous work 
involving both analytical calculations and simulations we 
and others analyzed the geometric effect of measuring 
the diffusion coefficient from the projected path of the 
protein [Mr- 26.. While these studies only included free 
diffusion, the lateral diffusion of an inclusion that inter- 
acts with the membrane shape is considered in recent 
studies [Ml |2"TH50"] . If the membrane shape fluctuations 
were not influenced by the protein, the effective diffusion 
coefficient would be increased compared to the free dif- 
fusion coefficient [231 1271 HH] • However, this simplifying 
assumption represents too severe an approximation. By 
including the backaction of the protein on the membrane 
fluctuations Naji et al. [30] showed that in equilibrium the 
effective lateral diffusion coefficient is reduced for which 
they were able to give an approximate expression. Parts 
of the current work are complementary to their study. 

On a more collective level, the interaction between 
membrane and embedded proteins can cause morpholog- 
ical changes. Leibler [31] studied a model with a protein 
density field that induces a spontaneous curvature capa- 
ble of causing an instability of the membrane. He finds 
two characteristic time scales for membrane fluctuations, 
with one of them potentially unstable in a certain wave 
vector range. Related work by Bivas and Meleard [32] 
on bending elasticity and fluctuations of spherical bilayer 
vesicles with additives reveals an additional characteris- 
tic time scales. In their work the bilayer membrane com- 
prises two individual sheets such that the third timescale 
results from the friction between the two layers. 

Divet et al. [33] studied an extension of Leibler's work 
that allows an exchange of proteins between the mem- 
brane and the surrounding fluid. Depending on the con- 
sidered length scale they find several relevant time scales 
for membrane height and density fluctuations. Surpris- 
ingly, this relevance of different time scales in membrane 
height correlations has, to the best of our knowledge, not 
been previously observed in experiments. 

In the present work we study the two interrelated ef- 
fects following from a protein-membrane interaction: the 
reduction of the lateral diffusion coefficient of the protein 
and the modifications of the static and time dependent 
height correlation functions of the membrane. In our 
study the additional energy caused by the insertion of 
the protein arises from an effective bending rigidity and 
spontaneous curvature of the protein. The dynamics of 
our system is dominated by two processes: the shape 
fluctuations of the membrane and the lateral diffusion of 
the protein. Taking into account the hydrodynamic in- 
teraction of a membrane with the surrounding fluid, we 
are able to derive a Langevin equation for the dynam- 
ics of the membrane shape, which becomes a function of 
the protein's position. Lateral diffusion of the protein is 
captured by another Langevin equation that takes into 



account the shape of and the interaction with the mem- 
brane. These two coupled equations of motions are the 
starting point for our analytical calculations and their 
numerical integration make out our simulation scheme. 

Using a systematic approach we present two main re- 
sults for lateral diffusion: we argue that in equilibrium 
the effective lateral diffusion coefficient of a protein that 
interacts with the membrane shape is universally de- 
creased compared to the free diffusion coefficient appli- 
cable if no interaction were present. Beyond this gen- 
eral argument we, furthermore, derive an explicit expres- 
sion for the effective diffusion coefficient of a protein with 
spontaneous curvature and bending rigidity by applying 
a path integral approach. To lowest order, our expres- 
sion agrees with that recently derived by Naji et al. [30] 
through an estimate of the power loss of the diffusing 
protein in the limit that the membrane shape minimizing 
the system's energy instantaneously tracks the protein's 
position. Our approach reveals that their expression re- 
sulting from a phenomenological approximation formally 
corresponds to the lowest order of an expansion in (/3k) _1 
using ft = (ksT)^ 1 with temperature T, Boltzmann's 
constant fee, and bending rigidity k of the membrane. 

We compare results of our simulation scheme with this 
lowest order expression and find good quantitative agree- 
ment. Several correlation functions contribute to the 
effective diffusion coefficient. The quantitative analysis 
of these contributions from our simulation results shows 
that the correlations between the response of the mem- 
brane and the preceding stochastic force acting on the 
protein effectively reduce the diffusion coefficient, while 
all other contributions would cause an increase. 

Concerning the altered membrane spectrum we analyt- 
ically develop an approximate expression for the height 
correlation function applicable for equal bending rigid- 
ity of protein and membrane. In the limit of slow pro- 
tein diffusion compared to membrane fluctuations we find 
two time regimes for the decay of height correlations: at 
small times the decay is dominated by membrane dynam- 
ics, while the diffusive time scale of the protein becomes 
the only relevant time scale at later times. Since we find 
this feature in the simulations not only for equal bend- 
ing rigidities of membrane and protein, we suggest that 
the experimental analysis of the late decay of dynami- 
cal membrane height correlations can provide a means to 
determine the effective diffusion coefficient of proteins in 
the membrane. To corroberate this argument we use our 
model to give a rough estimate suggesting that the effect 
should be visible in realistic systems. 

The paper is organised as follows: in the next sec- 
tion we introduce the model of the system and develop 
the equations of motion both for the protein and the 
membrane. In the limit of small ratios of temperature 
and bending rigidity these are then used in sees. |III A| 
and IHB to develop an exact analytical expression for the 
effective lateral diffusion coefficient and for the dynamical 
membrane height correlations, respectively. In sec. |IV A| 
we briefly explain our simulation scheme and motivate 
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the parameters used in the simulations. In sec. IV B we 



show that the derived analytical expression for the effec- 
tive diffusion coefficient shows good agreement with sim- 
ulations using parameters of typical experiments. We, 
furthermore, quantitatively analyse the correlation func- 
tions that contribute to the diffusion coefficient. The 
membrane height correlations are compared with our an- 
alytical expressions in sec. |IV C| In sec. [V] we show the 
determination of the effective diffusion coefficient from 
the late time decay of height correlations and discuss that 
this procedure should be experimentally feasible in real- 
istic system. We finally close with some conclusions. 



II. MODEL 

In our model we consider a single diffusing inclusion in 
a membrane with bending rigidity n that we describe in 
the Monge gauge. The small inclusion with radius a p has 
a spontaneous curvature C p and its stiffness may differ 
from that of the membrane by a factor of 7. The energy 
of the system of size L 2 may be expressed by 



H[h,K 

na p G(r 



2^ 



{(v^) 2 + 
)[ 7 (v^-c p ) 2 -(v^) 2 ]}, (1) 



with the height function h(r) that quantifies the dis- 
tance between the membrane and the position r on a 
flat reference plane. The particle position projected onto 
this plane is given by R. The function G(r — R) de- 
fines a weighting function of the particle that must ful- 
fil the normalization J L2 cPrG(r) — 1. In our sim- 
ulations we set the weighting function to a Gaussian 
G(r) — (na 2 ) -1 exp [— r 2 /a 2 ] such that the transition 
from membrane to particle is smooth. If we use the 
Fourier expansion h(r) = ^ k h(k) exp(zk ■ r) and 
h(k) = f L2 d 2 r h(r) exp(— ik-r) the Hamiltonian becomes 



H[Mk),R] = f ^5> 4 /i(k)M-k) 



2 L 2 



HT ft 

(7-1)^^ k 2 k' 2 G(-k - kV< k+k > R h(k)h(k')+ 



E fc 2 G(-kK k R Mk) \ + ^C 2 . (2) 



k,k' 

2 



L 2 



From this Hamiltonian it is in principle possible to nu- 
merically calculate the equilibrium height correlations 
(h(k)h(k')) applying methods used in [33]. While these 
methods are restricted to time independent equilibrium 
quantities, our simulation scheme, described later in the 
paper, allows us to not only obtain these quantities but 
also time dependent information. For the special case 
that both the protein and the membrane have the same 



bending rigidity, i.e. 7 
given by 

fc 4 



1, the height correlations are 



(h(k)h(-k)) 



1 



pira 2 p C 2 G(k)G(-k) 



(3) 



with the ratio of protein area to system size 

P = ™ 2 p /L 2 . (4) 

While eq. ([3| is derived for a single protein on the mem- 
brane the extension to several noninteracting proteins 
would lead to the same result with p resembling the 
overall area density of the proteins. Compared to the 
free membrane without protein, whose height correla- 
tions are given by the first term, the protein gives rise 
to an additive term that depends on the various parame- 
ters characterizing the particle. Following previous stud- 
ies [3TH33] . it is possible to define a fc-dependent effec- 
tive bending rigidity K e s(k) such that the spectrum of 
the membrane has the form of a protein-free membrane 
(h(k)h{-k)) bcc = L 2 /((3K cS (k)k 4 ). This leads to 

^1 = [l + p Kp na 2 p C 2 p G{k)G{-k)Y l (5) 

Note, that the addition of the protein in the membrane 
always leads to a reduction of the effective bending rigid- 
ity of the system as has been pointed out previously when 
inclusions are inserted into a membrane |31j . However, 
the effective bending rigidity cannot become negative, 
hence the membrane not unstable. 

Neglecting any geometric effects caused by the projec- 
tion of the protein path, which we have previously de- 
termined to be rather small for realistic membranes [24] , 
the diffusive motion of the protein and the thermal fluc- 
tuations of the membrane, i.e. the dynamics of the height 
modes h(k,t), are appropriately described by the follow- 
ing coupled Langevin equations: 

R(t) = -p P V R H + C(t) (6) 

with the stochastic forces £(t) and £(k, t) that are related 
to the mobilities p p = Dp/fceT of the protein [35] and 
A(fc) of the membrane, respectively, via the fluctuation- 
dissipation-theorems 

<G(*)> = 



and 



(6(*)Cm(tf)> = 2D p S hm S(t-i/), 



<£(M)> = 



(£(k,t)£(k',t')) = 2k B TA(k)L 2 Sv^S(t-t>). 



(8) 



(9) 



The mobility of the membrane takes into account the dy- 
namics of the membrane caused by the surrounding fluid. 
A hydrodynamical derivation involving the Oseen tensor 
leads to a mobility of A(k) = (Arjk)^ 1 for the undulations 
k ^ 0, with the viscosity r\ of the surrounding fluid [21]. 
For k = 0, i.e. the center of mass movement of the mem- 
brane, we set A(fe = 0) = since it does not influence 
the properties of interest in our study. 



4 



III. ANALYTICAL APPROACH 
A. Diffusion coefficient D e ff 

We derive an analytical expression for the effective 
diffusion coefficient D c g by exploiting that for biomem- 
branes typically (/3k) _1 <C f . 

We first determine the minimum of the energy ^ . The 



condition 



an 



dh(k) 



= leads to the equation 



= k 2 h_ k + 

(7- l)p^fc' 2 G(k + k')exp[i(k + k') .B]hv 

k' 

+ 7 C>a 2 G(k) exp(ik • R) , (10) 

for the height modes that minimize the energy. Using 
the ansatz — -pr exp(— ik-R) the energy is minimal 
for 



Bk = -7<W£M J ^G(q) , 



with the matrix 



M k , q = 6 Kci + (7 - l)p G(q + k) 



(11) 



(12) 



Inserting this result into the Hamiltonian shows that the 
energy minimum does not depend on the particle position 
as expected from the isotropy of the particle position. 

The first question we will address is whether the ef- 
fective diffusion constant is larger or smaller than the 
free diffusion coefficient D p applicable without coupling, 
i.e. 7 = 0. The degrees of freedom in the Hamiltonian |2| 
are given by the membrane modes h(k) and the position 
R of the protein. Due to the appearance of the mem- 
brane modes up to quadratic order in the Hamiltonian, 
the thermal averages (|/i(k)| 2 ) remain bounded. The 
position of the protein, however, is not bounded, such 
that diffusive motion is possible. Effectively, the protein 
moves in a time-dependent periodic potential given by 
the height modes h(k,t). Diffusion in periodic potentials 
has been previously considered in a large number of stud- 
ies. If the particle is only subject to the potential and no 
other external force it is easily shown, that the effective 
diffusion coefficient of the particle is always smaller than 
or equal to the free diffusion coefficient [5S1 [37] . Thus 
we conclude that the effective diffusion coefficient D c s of 
a protein whose interaction with the membrane depends 
on the shape obeys 



D cS < £> p 



(13) 



in all situations without external driving forces. 

While we know that the diffusion coefficient is in gen- 
eral reduced due to the membrane we will derive an ex- 
plicit analytical expression to quantify the effect for our 
model, the validity of which we will discuss and analyze 



through simulations. A quick glance at the Langevin 
equations ([6| and Q shows that they are highly nonlin- 
ear such that the exact solution is not straightforward. 
In order to develop an expression for D c s we must, there- 
fore, apply certain approximations. 

In a path integral description [38] the probability dis- 
tribution V of the paths R(f ) of the diffusing protein and 
of the height modes h(k, t) follow from the weight of noise 
fluctuations and are given by the functional 



V[K(t),h(k,t)] ~ exp 
with the function 



~ J*drL(R(T),h(k,T)) 



(14) 



L(R,h(k)) = (r + /! VrH) 2 



E 



i 



2A(fc) 



h(k) + A(k) 



dh(k) 



If we introduce the deviation 



y* = h(k) - h k (R(t)) , 



(15) 



(16) 



of the membrane shape from the instantaneous mem- 
brane shape that minimizes the energy of the system the 
Hamiltonian may approximately be written in the form 



1~L — 



S 2 H 



2/k' 



(17) 



MR) 



where T-Lq = H[h(t)] is the energy minimum. Since the 
second functional derivative of the energy with respect 
to the height is proportional to the bending rigidity n, it 
is convenient to define 



6 2 H 



6h(k)Sh(k') 



(18) 



MR-) 



Up to total derivatives leading only to boundary terms 
we can then rewrite the function L as a function of R 
and j/ k 



L(R,y(k)) 

y— 

^ 2A(fc) 



1 

2^o ' 



R 



(R-VrK +2y k (R-V R )/i k + |y k | 



Mo(/3k) 2 ( ^(ykV R U k ,k'2/k' - 2/kUk,k'V R /ik') 



k,k' 



i(/^) 2 £A(k)(]TUk,k' 2 /k< 



(19) 



To determine the effective mobility of the particle it 
would be necessary to integrate out the deviation y in the 
probability distribution V . Since this cannot be done ex- 
plicitly, we employ a saddle point approximation, i.e. we 
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replace all the possible paths y k (i) with the path yu(t) 
that minimizes the function L and hence contributes to 
the probability distribution (14 1 the most. The path 



yw(t) follows from the Euler-Lagrange-equations 



d dL 
dt dy k 



dL 



(20) 



As before, we assume that only small variations y k in the 
height are relevant allowing us to linearize the resulting 
differential equations 

j^jUt) - (w 2 E A ( k ")(E^"^'^)^.^ 



k",k' 



1 d 2 /l k (R(i)) 

A(k) dt 2 



(21) 



or, in a simplified notation, 



m - cwe^s* = , (22) 

k,k' 

using the positive definite matrix A kk /. The homoge- 
neous solution of this equation is a simple relaxation on 
time scales proportional to {/3k) _1 and plays, therefore, 
a minor role for large values of (3k. We are now inter- 
ested in the significance of the bending rigidity k on the 
inhomogeneous solution. To this end we drop all depen- 
dencies of k in eq. ( 22 ) as though the system only had a 
single wave mode without altering the order of (3k. The 
inhomogeneous solution is then given by 



yinh(i) 



2(3kVA 



dr 



dT 2 ' 

For large (3k and slow protein diffusion the mem- 
brane shape minimizing the energy h(R(t)) only weakly 
changes on the relaxation time scale. Hence, membrane 
shape deviations y k are of the order 0(((3k)~ 2 ). Inserting 
this result into eq. ( 19 ) and keeping only leading orders 



of ((3k) we arrive at 
L(R,S(k)) = 

^ 2 + E^*-Va 



2fiQ 



0((/3k)- 



(24) 



It is now possible to identify an effective diffusion coeffi- 
cient -D e ff for the diffusing protein from the prefactor of 
the R 2 term 



Dq_ 



= i+E 

k 

= i+E 



Mo 



V R /i k 



Mo ^k^-k 

A(fc) k 2 



(25) 



with B k from eq. (11) and using isotropy in the x- and 



y-direction of the system. This systematic derivation 
of an analytical expression for the diffusion of a protein 
that interacts with the membrane shape constitutes our 
first main result. The analysis of experiments with both 
model and biological membranes, shows, that (/3k) -2 is 
typically smaller than 0.01, hence, sufficiently small to 
expect a wide applicability of this expression. The first 
line of this expression agrees with the result of Naji et 
al. [30]. In their derivation they apply an adiabatic ap- 
proximation assuming the membrane shape to instanta- 
neously follow the path of the protein. The effective dif- 
fusion coefficient is then derived from an estimation of 
the power loss of the diffusing particle. Our derivation 
identifies their approximate result as the lowest order of 
an expansion in (/3k) _1 . 

After the general solution for the effective diffusion co- 
efficient we will now turn to the special case of a pro- 
tein with a weighting function expressed through Dirac's 
delta- function G(R — r) = 5(R — r). In this case the 
Fourier-transform of G is independent of the wave- vector 
k, such that G(k) = 1. This leads to the height-mode 

fe k (R)= 1 p expHk-R), (26) 

k 2 1 + (7- IJpEk' 1 

minimizing the free energy %. The resulting energy min- 
imum is 

H[h} = ^na 2 p C% — 1 . (27) 

Following the above procedure but now inserting the spe- 
cial choice for the weighting function leads to the effective 
diffusion coefficient 



D 



Mo 1 



7 2 C> 2 < 



- + ^A(fc)fc 2 (i + ( 7 -i WEk ,i)) 2 ' - 



In a next step we will evaluate the sums over all possible 
wave- vectors k. On the one hand the smallest value of the 
x or y component ist defined by the system size, k x>la i a = 
2n/L, and is thus approximately zero for large systems. 
On the other hand the largest value fc Xjmax = 2tt/cl c is 
limited by a microscopic cutoff length a c that corresponds 
to the size of the lipids. Thus the expression pX) k 1 m 



the denominator of eq. (28) leads to 7ra 2 /a 2 using the 
definition of p given in eq. Q. The evaluation of the 
other sum over k depends on the specific form of the 
membrane mobility. Using A(k) = (Arjk)^ 1 makes the 
evaluation of ^ k fc _1 necessary, such that the expression 
for the ratio of the diffusion coefficients becomes 



= 1 



^ 47?7 2 C 2 7ra^ 
1 + (7-l)| 



In 



V2 + 1 \ L 2 
V2-I) a c 



(29) 



Typically the stiffness of a diffusing protein will be sig- 
nificantly larger than that of the membrane. In the limit 
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7> 1 the resulting effective diffusion coefficient is given 
by the relation 



Do 



1 + HoA-qClTT In 



2„3 



V2- 1 



(30) 



Equation ( 28 ) shows that a non-zero spontaneous cur- 



vature C p is crucial in order to have an influence on the 
diffusion coefficient of the protein. Within our model a 
mere difference in the bending rigidity of the membrane 
and the protein, i.e. 7^1, does not lead to an altered 
diffusion coefficient. 

Before testing our expression by comparing it with sim- 
ulations, we will discuss two limiting cases for eq. (28 1. 



Without altering the general conclusions we will give 
this discussion only for a single fc-mode. As was ex- 
pected, our expression reveals an effective diffusion co- 
efficient that always has an upper bound of D . For very 
small ratios /io/A (A = A(k)), i.e. if the membrane is 
much more mobile than the protein, the reduction of 
the effective diffusion coefficient is linear in this ratio 
D e g/D « 1 — |£?k| 2 Mo/A. The free mobility of the pro- 
tein dominates its effective movement. The influence of 
the membrane is weak since it can adjust quickly to the 
position of the protein. In this situation our approxima- 
tion that the system's energy is always close the mini- 
mum is fulfilled and our expression will serve as a very 
good estimate. For the limit in which the membrane 
moves much slower than the particle, which corresponds 
to the scenario of diffusion in an (almost) fixed periodic 
potential, our expression predicts the asymptotic behav- 
ior D e g/Do |£?k|~ 2 A//io- Thus the movement of the 
protein is mainly dominated by the membrane mobility 
such that it is strongly slowed down. However, if the dif- 
fusing particle effectively sees a fixed energy landscape 
the stochastic motion enables the protein to hop from 
one energy minimum to another, such that our previous 
approximation that the protein always stays very close to 
the position of the instantaneous energy minimum may 
no longer be valid. 



Langevin equation ^ is then given by 
h{k,t) =h(k,0)e' t/TM{k) +e~^™^ 



dt 1 e*' /™ {h) 



k 2 T M {k) 

(31) 



with the fc-dependent membrane time scale 

TMi^^^/iKk 3 ). (32) 

Using the fluctuation-dissipation theorem ^ and as- 
suming that the particle diffuses on timescales much 
larger than rjv/(fc) the height correlation function is given 
by 



(MM)M-M')) = 



L 2 



e -\t-t'\/T M (k) + 



@Kk 4 

PK P ira 2 p C p 'G(k)G{--k)e- lt - t ' l/TD{k) 



with the diffusive time scale 



TD (k) = (D cS k 2 y 



(33) 



(34) 



In order to arrive at eq. (33) we use (exp[— ik • (R(<) — 
R (^))]> = e - l t-t 'l/ TD W which follows for diffusive mo- 
tion of the protein with an effective diffusion coefficient 
D cB . 



Equation (33) shows that the dynamics of the height 
correlation function of the membrane is determined by 
the two timescales present in the system: the membrane 
time scale Tjvf(fc) and the diffusive time scale Tjj(k). Our 
calculations assumed tm(&) < Tu(k), thus, the decay of 
height correlations for small times t—t' will be dominated 
by the membrane dynamics while for large times the dif- 
fusion of the particle takes over. While our calculation is 
strictly valid only in the case of 7 = 1 the qualitative be- 
haviour persists also for the case 7 7^ 1 as will be shown 
when we present simulation results. Note, that the naive 
usage of an effective binding rigidity K c e{k) would lead to 
a single timescale for each mode k. Since the properties 
of the system are clearly dominated by two timescales 
the concept of K e fi{k) is only applicable for properties 
that are not time dependent. 



B. Temporal decay of membrane height 
correlations 



IV. SIMULATIONS 



Using the equations of motion ([6| and Q , a calculation 
of the full height correlation function (/i(k, t)h(k! , t')) is 
not feasible analytically. We rather determine this quan- 
tity from our simulation scheme. Nevertheless, in order 
to gain an understanding of the possible contributions 
to the correlation function, it is instructive to consider 
the special case of equal bending rigidities of the particle 
and the membrane, 7 = 1. The general solution of the 



A. Scheme 

Our simulation scheme comprises the numerical inte- 
gration of the two coupled Langevin equations (6]) and 
([7]). However, the equation of motion given for the pro- 
tein in eq. §6§ neglects that the particle actually dif- 
fuses along the membrane, in other words a curved sur- 
face. The shape of the membrane influences the Langevin 
equation, the exact form of which is given in refs. [231 12"7] 
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and used in our simulations. Thus the free diffusion co- 
efficient Dq used in the simulations is slightly larger than 
the value of D p in cq. (|6| [55] . 

The membrane is mapped on a square N x AT-lattice 
such that the length of the system is L = N £ with the 
lattice spacing t. The membrane shape is evolved in time 
by a time discrete version of eq. ^ in Fourier space. This 
part of the scheme is an extension of the Fourier Space 
Brownian Dynamics simulation method introduced by 
Lin and Brown [MEET] . 

After every update of the membrane shape the posi- 
tion of the particle is altered by using a discrete version 
of eq. (|6j). However, the particle's position is not evolved 
on the lattice. The membrane height at the particle posi- 
tion that enters in the equation of motion is determined 
through linear extrapolation from the height at the four 
nearest neighbour lattice sites. The shape of the mem- 
brane in real space h(r, t) is determined by use of fast 
Fourier transforms implementing the FFTW library |42j . 
For a more detailed account of the simulation scheme we 
refer the reader to refs. [55J [27] . 

All simulation results presented in this paper were per- 
formed on a 64 x 64-lattice with a lattice spacing of I = 
lOnm. The radius of the protein is set to a p = 21. The 
fluid surrounding the membrane is water with a viscosity 
of r? = lfr 3 kg/(ms) or rj = 2.47 x 1CT 7 's/ '(fit 2 ) in the 
units of our model at T = 300!^. The discrete integra- 
tions of both membrane shape and particle position are 
performed with a timestep of At = 10 _9 s that is signifi- 
cantly smaller than the smallest time scale TM,xam in the 
system. If not stated otherwise the diffusion coefficient of 
the protein is set to D = 5 x 10~ 8 cm~ 2 /s= 5 x 10 4 ^ 2 /s. 
This ensures that the membrane time scale is always 
smaller than the diffusive timescale tm(&) < Td(^) as 
is the case in real biological systems. Simulation runs 
were performed with 8 x 10 6 integration steps resulting 
in trajectories that last for 8ms, which is approximately 
five times the longest membrane time scale TM,max- 

The 

graphs presented in the following are the results of av- 
eraging over a minimum of 500 independent trajectories, 
where the first 10 6 timesteps were not taken into account 
in order to ensure equilibration of the membrane configu- 
ration and the particle position relative to the membrane 
shape. 



B. The effective diffusion coefficient 



To test our explicit expression ( 25 ) we have performed 



elaborate simulations using the scheme described in the 
previous section. In fig. [I] we present the resulting 
D e ff/Do as a function of pj with the ratio of protein area 
to system size p, eq. Q , and the ratio of the protein to 
membrane bending rigidity 7 for three different protein 
mobilities. The detailed parameters of the simulations 
are given in the figure caption. The comparison of the 
simulation results and the corresponding analytical ex- 
pression shows very good agreement for all the chosen 
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FIG . 1 : Ratio D e g / Do as a function of the coupling coefficient 
P7 for the given free diffusion coefficients Do in units of £ 2 /s. 
The bending rigidity of the membrane is /3k = 5 and the 
spontaneous curvature of the protein C p t = 1. 
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FIG. 2: Force correlations integrated twice in time to show 
the respective contribution to the mean squared displacement 
of the protein for 7 = 1 and C v l = 1. The results apply for 

T $ T . 



parameters. Since the simulations were all performed 
with (/3k) -2 = 0.04 we are well within the limits, where 
we expect our analytical result to hold. These simulation 
parameters were chosen because they represent realistic 
parameters for biological systems. We thus conclude that 
our approach to determining the reduction of the effec- 
tive diffusion coefficient will be of use in experimental 
studies. 

In the following we will quantitatively analyse the con- 
tributions of the correlation functions entering the effec- 
tive diffusion coefficient. The mean squared displacement 
(AR 2 (<)) for the diffusing protein is formally given by in- 
tegrating eq. ^ twice in time 

(AR 2 (t)) = f dr f At' <9 t R(t) • 9 T 'R(r')> 
Jo Jo 

= J* dr J* dr' [ M 2 (f (r) • f (r')) + (C(r) ■ C(r')> 

+Mp<f(r).C(rO}+Mp<C(r)-f(rO>] , (35) 
with the conservative force f(i) = — VrK[/i, R] that the 



membrane exerts on the protein. Thus the mean squared 
displacement has several additive contributions. Since 
the effective diffusion coefficient follows from the slope of 
the mean square displacement as a function of time via 
(AR 2 ) = 4D^ SD i, also the diffusion coefficient has vari- 
ous additive contributions. In fig. [2j we display the sim- 
ulation results for the various parts of the mean squared 
displacement for a chosen set of parameters. The corre- 
lations of the stochastic force acting on the particle lead 
to a ADpt behavior as is expected from the fluctuation 
dissipation theorem of eq. (|8j). As we have argued before, 
the slope of the particle's mean squared displacement as 
a function of time is smaller than 4D p , hence, one of 
the additive terms must be negative. However, the force 
correlations (f(r) • f(r')) obviously also lead to an addi- 
tive contribution, which we find to be quite small for the 
parameters of our simulations. Due to causality correla- 
tions (f(r) • C( r ')) w hh r ^ t 1 must be zero, such that 
the remaining correlations (C(t') • f(V)) are the cause of 
the reduction of the diffusion coefficient, as we clearly see 
from the simulation results. This contribution expresses 
the reaction of the membrane to the random force acting 
on the protein: If the random force moves the particle 
during a small discrete timestep, the interaction of the 
protein with the membrane will slightly change the shape 
of the membrane during the next timestep such that the 
system comes closer to the energy minimum. However, 
it cannot be reached during such a short time, leading 
to the membrane "pulling back" the protein to its initial 
position before the random movement. This explains the 
sign of the corresponding correlation function. An impor- 
tant aspect here is that the membrane shape reacts to the 
movement of the protein. If the membrane shape evolves 
independently from the particle position these correla- 
tions do not exist leading to an increase in the effective 
diffusion coefficient [241127] , 
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FIG. 3: Mean squared height of a membrane with protein 
normalised by the height correlations of an unperturbed mem- 
brane as a function of k for equal bending rigidity of the pro- 
tein and the membrane, 7=1. The values for the membrane 
bending rigidity /3k and the spontaneous curvature C v are 
displayed in the legend. We define E(k) = (j3nk 4 ')~ 1 . 
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FIG. 4: Effective bending rigidity « e ff(fc) as a function of k for 
the given values of 7. The bending rigidity of the membrane 
is /3 k = 5 and spontaneous curvature of the protein C p £ = 1. 



C. Membrane height correlations 

In the following we will elucidate that our simulation 
results for equal bending rigidity of membrane and pro- 
tein, 7 = 1, agree very well with the analytical expres- 
sions given in eqs. ([3| and (33 1. Furthermore, we will 
show that the qualitative features of these equations are 
also observed in the more general case 7 7^ 1. 

In fig. [3j we present the height correlation spectrum 
(h(k)h(— k)} as a function of k 4 for 7 = I and different 
membrane bending rigidities k and spontaneous curva- 
tures C p of the protein. In order to focus on the influ- 
ence of the protein we have normalised (h(k.)h(— k)) by 
the spectrum of a bare, protein-free membrane. While 
the symbols represent results from the simulations the 
solid lines follow from eq. ([3| using the Gaussian weight- 
ing function G(r— R) given in sec. [TlJ The influence of the 
protein is most pronounced for small wave vectors k or 
large length scales and decreases with increasing k to the 
value of the bare membrane without protein. Membrane 



fluctuations on length scales significantly smaller than 
the inclusion's size are not influenced by the interaction 
of the protein with the membrane. Comparing simula- 
tion results with the analytical expression ^ we find that 
the agreement is very good as was of course expected. In 
fig. [4j we plot the effective bending rigidity n c s{k) as a 
function of k as determined from the height correlations 
that result from simulations with a constant spontaneous 
curvature C p and membrane bending rigidity k, but dif- 
ferent bending rigidity ratios 7. For 7 = 1 the simula- 
tion result agrees very well with eq. [5j With increasing 
the stiffness of the particle we find that the qualitative 
behaviour remains similar, however, n e g is even more re- 
duced. Effectively, an increase in the bending rigidity of 
the protein leads to a softening of the system. Our results 
indicate that K c g(k) saturates with increasing 7. To cor- 
roborate this assumption simulations with even higher 7 
would need to be performed, but have turned out to be 
very demanding computationally. 

We now turn to the temporal decay of height correla- 
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FIG. 5: Height correlation functions {h(k, t)h(— k, 0)} for two 
arbitrarily chosen fc-values as a function of time t. The slopes 
of the solid lines determine the two time regimes dominated 
by TM(fc) for small t and ro(fe) for large t. 



tions. In fig. [5] we display (/i(k, t)h(— k, 0)) as a function 
of time for two arbitrary wave numbers k. The main 
plot considers 7 = 1, the inset 7 = 6. Equation (l33|) 



suggests that the two relevant time scales in the system, 
that are well separated in our calculations, become ob- 
servable: for small times the decay is dominated by the 
membrane time scale tm(^), while for larger times the de- 
cay is predominantly influenced by the movement of the 
protein and hence the corresponding time scale is ri)(fc). 
Regarding the simulation results we see indeed that the 
behaviour of the correlations is dominated by a fast de- 
cay at small times and a slower decrease for large times. 
While eq. ( [33] ) is an approximate result only for 7 = 1 
we find that this feature of two dominating time scales is 
qualitatively also observed for 7^1. The quantitative 
fit of the small time behaviour with an exponential de- 
cay with the characteristic time t m(^) shows very good 
agreement for both considered values of 7. For 7 1 the 
Hamiltonian of the system ^ causes an additional con- 
tribution to the inverse characteristic time that depends 
on 7. However, for the parameters of our simulations this 
contribution is obviously negligible. 



At large times the decay is expressed through the time 
scale Toik) that is a function of the effective diffusion 
coefficient D c g of the protein along the membrane. The 
typical method to identify D e f{ is to regard the temporal 
evolution of the mean squared displacement of the pro- 
tein using the relation (AR 2 ) = 4Z?^ SD i, as explaned 
above. Using the so determined value of the diffusion 
coefficient we find for large times that the results in 
fig. [5] are well approximated by an exponential behaviour 
~ exp[— t/TD(k)} for both 7 = 1 and 7 = 6. Overall, 
for 7 = 1 we observe that the behaviour of the height 
correlations is well described by eq. (33). 



FIG. 6: Diffusion dominated decay time Tr>(k) as a function 
of k for 7 = 1 and the given values for j3n and C v . Sym- 
bols represent simulation results derived from fitting height 
correlation functions at later times; solid lines display the 
theoretical time scale using the effective diffusion coefficient 
determined through the mean squared displacement of the 
protein. In the inset the effective diffusion coefficient deter- 
mined from the long time decay using D^j = (fc 2 TD(fc))~ 1 is 
plotted as a function of k (symbols) . The horizontal lines give 
-Deg^ D determined from the mean squared displacement. 



V. DIFFUSION COEFFICIENT EXTRACTED 
FROM MEMBRANE SPECTRUM 

A. Determination of D e fr from simulations 

In the discussion of fig. [5] we used the mean squared 
displacement to determine the diffusion coefficient of the 
protein. However, the exponential decay of large time 
height correlations offers an alternative method to ex- 
tract D c ff. The Tr>(k) resulting from exponential fits 
to the late time decay as a function of k are plotted 
in the main graph of fig. [6] for 7 = 1 but different 
k and C p . Using the previously determined value for 
the diffusion coefficient we find a good agreement with 
T D (k) = (/c 2 D^ SD ) _1 (solid lines). Thus without prior 
information on the mean squared displacement of the 
proteins, it is possible to identify D c s from the height cor- 
relations of the membrane using the behaviour of To{k) 
as a function of the wave- number k. In the inset of fig. [6] 
we plot D^g = (fc 2 T£)(fc))~ 1 as a function of k. We find 
that the resulting diffusion constant agrees very well with 
the diffusion constant Z?^ SD . 

While fig. [6] only considers 7 = 1, we will now show 
that the characteristic timescales tm and tjj can also be 
identified for 7 7^ 1. The inset of fig. [7] displays rjvf(fe) 
as a function of k as determined from the exponential 
decay of the height correlations at short times for differ- 
ent values of 7. Apart from 7 the other parameters of 
the membrane and the particle are kept constant. We 
find that the results do not depend on the rigidity of the 
protein. The dominant time scale for short times is only 
determined by the properties of the membrane and co- 
incides with the correlation time of a freely fluctuating 
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FIG. 7: Diffusion time rn(k) as a function of k for the given 
values of 7. All results apply for C v l = 1 and /3k = 5. The 
symbols result from simulations while the solid lines are given 
by (_D^ SD fc 2 ) -1 . The inset shows the corresponding mem- 
brane time scale tm (k) determined by fitting the initial decay 
of the height correlation function (symbols) and given theo- 
retically (solid line). 
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FIG. 8: Crossover time r co (fc) from membrane fluctuation to 
diffusion dominated decay of height correlations as a function 
of wavenumber k, see eq. p6|, for 7 = 1, D c fi = 10" 8 cm 2 /s, 



y/na p C p = 1, and /3k — 10. The dotted line corresponds to 
the typical resolution time of video microscopy. 



membrane without protein. In the main plot the late 
time diffusive time scale td is plotted as a function of 
k. These results clearly depend on the rigidity of the 
protein, but only because the diffusive motion is influ- 
enced by the protein-membrane interaction. If we de- 
termine the effective diffusion coefficient from the mean 
square displacement and compare the thus calculated tjj 
(solid lines) with the simulation results from the late time 
exponential fits to the height correlations (symbols) the 
agreement is again very good. Thus, if we had not had 
the possibility to determine the mean squared displace- 
ment of the protein we could have determined D c g solely 
from the time dependence of the height correlations. 



than approximately 4pm . 

If we have a general time-dependent function of 
the form Aexp(-at) + B exp(-bt) with a ^> b a 
good estimate for the crossover time is given by (A + 
B) exp(— ar co ) = B exp(— 6r co ). Within our model the 
crossover time for 7 = 1 is given by 



m(fc) r D (k) 



x In 



(f3 K pira 2 p C*G(k)G(~k)y 



(36) 



B. Estimate of experimental feasability 

In the following we will show that our suggested 
method to determine the lateral diffusion coefficient of 
proteins from height fluctuations in a membrane is ex- 
perimentally feasible. Height correlation functions can 
be determined by video microscopy as explained in the 
introduction. If T co (fc) is the crossover time from the 
decay of correlations caused by the membrane dynam- 
ics to that dominated by the diffusive timescale of the 
proteins, experiments must meet two conditions in or- 
der for the crossover to become observable: on the one 
hand T co (k) must be larger than the temporal resolution 
of the camera used in the experiment. Values given in 
previous studies [15] are on the order of 0.03s. On the 
other hand r co (k) must be smaller than the experimen- 
tally accessible total time that is on the order of minutes. 
Since the crossover time is a function of the wave vector 
the restrictions for r co (k) must, furthermore, be fulfilled 
for experimentally accessible wave numbers. The spatial 
resolution of video microscopy allows for k- values smaller 



While lateral diffusion coefficients of proteins in mem- 
branes are on the order of 10~ 8 cm 2 /s [IT], the spon- 
taneous curvature is not so well determined. In fig. [8j 
we display the crossover time T co (k) as a function of 
the wave number for different protein densities p and 
D = 10~ 8 cm 2 /s, (3k = 10, and */ira p C p = 1. For the 
regarded protein densities we find that for small wave 
numbers t co becomes larger than the typical temporal 
resolution of experiments. Wc find that for wave numbers 
that lie within the experimental range it should be pos- 
sible to observe the two characteristic time scales. Thus, 
the determination of the lateral diffusion coefficient from 
the diffusion dominated decay of membrane height corre- 
lations should be feasible. However, the interesting wave 
number range is reasonably narrow, since t co is strongly 
increasing for smaller k. Note, that eq. (36) is only valid 



for the situation of equal bending rigidity of the protein 
and the membrane. In general this is obviously not the 
case, however, the height correlations displayed in fig. [5] 
for 7 ^ 1 let us assume that the crossover time is only 
weakly influenced by the bending rigidity of the protein 
and that our estimate remains valid. 
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VI. CONCLUSIONS 

In this paper we have considered the influence of a 
protein interacting with a fluctuating membrane via its 
bending rigidity and spontaneous curvature on both the 
dynamics of the protein and the membrane. The quan- 
tities we have looked at in detail are the lateral diffusion 
coefficient of the protein and the height correlations of 
membrane fluctuations by use of analytical calculations 
and Langevin simulations. We argue that the lateral 
diffusion coefficient of a protein that interacts with the 
membrane is always reduced compared to its bare diffu- 
sion coefficient as long as there are no external driving 
forces or active processes. Using a path integral approach 
we could derive an analytical expression for this reduc- 
tion that is valid within the lowest order of a 
cxpansion. Our simulations with parameters that resem- 
ble those of real experiments show a wide applicability 
of this expression. In addition a closer look at the cor- 
relation functions that contribute to the reduction of the 
diffusion coefficient shows that the correlations between 
the stochastic force acting on the protein and the re- 
sponse of the membrane to the movement of the protein 
are responsible for the reduction. 

The diffusion of the protein is obviously correlated with 
the height correlations of the membrane. The determina- 
tion of the height correlations for the case of equal bend- 
ing rigidity of the protein and the membrane reveals that 
the protein-membrane-interaction has a significant influ- 
ence compared to a free membrane. The most predom- 
inant feature is that the temporal decay of correlations 
does not only display the timescale one would expect 
from the membrane, but that the diffusive time scale of 



the influencing protein becomes important. In realistic 
biomembrane systems these two time scales are well sep- 
arated such that a crossover from the initially fast decay 
of membrane fluctuations to the slower protein diffusion 
dominated decay, is likely to be observed in experiments. 
Since the decay at later times is directly related to the ef- 
fective diffusion coefficient of the protein, we suggest that 
the measurement of membrane fluctuations might actu- 
ally provide a means to determine the lateral diffusion 
coefficient of the inserted proteins. 

Our systematic approach to lateral diffusion of a pro- 
tein interacting with the shape of the membrane and the 
related influence on the membrane fluctuation spectrum 
can be extended in various directions. The first ques- 
tion arising from our study is to work out corrections 
to the (/3«;) _1 -expansion and to estimate their relevance. 
A further perspective resulting from our analysis is the 
interesting limit A//io — > when the protein effectively 
moves in a fixed membrane configuration. This situation 
is interesting theoretically, since the diffusing particle no 
longer "drags along" the membrane, but is hindered in 
its movement by potential barriers caused by the inter- 
action of the particle with the membrane. Finally, while 
we have so far only considered membranes that are on 
average flat, the extension to ruffled membranes poses 
an interesting challenge with significant relevance for bi- 
ological membranes like the endoplasmic reticulum or the 
cristae in mitochondria. 
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